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FOREWORD 


This  publication  was  prepared  under  Task  9,  “Forecasting  Rules  for  Specified  Strategic 
and  Tactical  Areas,”  and  describes  a  method  for  obtaining  the  maximum  or  limiting  probabil¬ 
ities  for  predictions  made  by  standard  quantitative  techniques.  The  method  has  a  theoretical 
basis,  but  this  will  not  be  emphasized  here.  Instead,  the  calculations  of  the  quantities  presented 
in  this  report  are  explained  and  interpreted,  with  examples  of  how  they  can  be  applied. 

The  limiting  probabilities  are  concluded  to  be  valuable  in  determining  the  reliability  of  the 
results  produced  by  the  given  forecast  techniques. 

This  report  was  written  by  Dr.  Thomas  A.  Gleeson,  Navy  Weather  Research  Facility  con¬ 
sultant  from  Florida  State  University.  Liaison  with  Dr.  Gleeson  during  the  preparation  of  this 
report  was  done  by  the  Task  Leader,  Mr.  S.  Donald  Case.  Mr.  Robert  S.  Haltiner  performed 
the  mathematical  review  and  Mr.  John  M.  Mercer  the  final  edit  of  this  report  for  the  Navy 
Weather  Research  Facility, 

This  publication  has  been  reviewed  and  approved  on  10  March  1966  by  the  undersigned. 


Commander,  U.S,  Navy 
Officer  in  Charge 

U.S.  Navy  Weather  Research  Facility 
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1.  INTRODUCTION 


In  any  large-scale  meteorological  situation, 
there  is  no  more  than  a  finite  number  of  synoptic 
measurements  available  for  map  analyses.  The 
lack  of  information  between  these  reports  makes 
such  analyses  somewhat  uncertain  and  implies 
that  forecasts  based  upon  them  will  also  be  un¬ 
certain,  to  some  extent  at  least.  As  one  meteor¬ 
ologist  [8]  puts  it,  “A  forecast  of  weather  can 
never  be  more  than  a  statement  of  probabilities.*  * 

A  second  consequence  of  the  inevitable  data 
void  between  reports  is  that  because  small  fea¬ 
tures,  which  may  be  significant  for  predictions, 
are  smoothed  out  of  the  analyses,  any  evaluation 
of  the  number  or  sizes  of  analysis  errors  will 
tend  to  be  underestimated.  In  other  words,  we 
can  do  no  better  than  to  minimize  the  errors  and 
thereby  maximize  the  estimated  probabilities  of 
success  of  large  scale  predictions.  We  conclude 
that  probabilities  are  necessary  ,  and  maximum 
probabilities  are  sufficient  in  synoptic  analysis 
and  forecasting. 

The  purpose  of  this  publication  is  to  des¬ 
cribe  a  method  for  obtaining  such  maximum  or 


limiting  probabilities  for  predictions  made  by 
standard  quantitative  techniques.  The  method 
has  a  theoretical  basis,  but  this  will  not  be  em¬ 
phasized  here.  In  application  to  standard  fore¬ 
cast  formulas  (extrapolation,  barotropic  models, 
etc.),  partial  derivatives  must  be  taken  to  arrive 
at  equations  needed  for  probability  predictions. 
But  when  the  latter  equations  are  at  hand,  their 
practical  use  requires  nothing  more  than  alge¬ 
bra. 

In  what  follows,  no  consideration  is  given 
to  errors  in  synoptic  reports  themselves,  nor 
to  imperfections  of  standard  forecast  techniques. 
This  is  another  reason  why  the  computed  prob¬ 
abilities  of  success  become  maximized. 

In  presentation  of  the  method,  it  is  deemed 
desirable  to  use  frequent  examples,  beginning 
with  the  simplest.  This  is  done  in  the  next  sec¬ 
tion.  A  summary  of  the  method  is  given  in  chap¬ 
ter  3,  and  subsequent  sections  are  devoted  to 
discussion  and  other  applications  of  the  general 
procedure.  References  [2,  3,  4,  5,  and  6]  can  be 
consulted  for  further  details. 
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2.  EXAMPLES  OF  METHOD 


2.1  Locotion  of  an  Isobar 

Because  of  the  absence  of  information  be¬ 
tween  synoptic  observations,  the  analysis  of  a 
field,  sea-level  pressure  for  example,  will  be 
in  error  to  some  extent.  Figure  2.1  shows  an 
isobaric  analysis  from  a  portion  of  a  hypotheti¬ 
cal  pressure  field.  The  1010-mb.  isobar  is 
drawn  through  an  arbitrary  point,  A  ,  in  this  case. 
However,  unless  there  happened  to  be  a  station 
at  A  reporting  1010  mb.,  the  true  pressure  there 
probably  would  not  be  exactly  this  value,  indi¬ 
cating  the  isobar  to  be  misplaced. 


Figure  2.1.  Example  of  Isobar  Pattern. 


If  we  define  an  x  axis  through  A  and  perpen¬ 
dicular  to  the  isobars,  as  infigure  2.1,  the  error 
in  location  of  the  1010-mb.  isobar  can  be  repre¬ 
sented  schematically  by  the  interval  Ax' .  Of 
course,  the  magnitude  of  Ax  cannot  be  found  from 
the  analysis  because  continuous  observations  are 
not  available.  But  it  is  possible  to  say  something 
about  the  probability  that  a  specified  magnitude 
of  Ax  will  exist. 

We  can  regard  the  synoptic  reports  as  being 
located  at  random  relative  to  the  true  pattern 
prevailing  at  a  given  time  in  a  given  geographic 
region.  (This  is  not  the  same  as  saying  the  re¬ 
ports  are  at  random,  geographically.) 


the  isobar  at  point  A  .  Now,  suppose  the  x  axis 
to  be  fixed  relative  to  the  true  pattern,  and  we 
consider  other  locations  of  the  1010-mb.  isobar 
along  x ,  as  would  be  obtained  from  independent 
analyses  of  the  same  basic  pattern  but  for  differ¬ 
ent  random  arrangements  of  synoptic  reports. 
In  all  such  arrangements,  assume  that  the  aver¬ 
age  density  of  the  observing  networks  remains 
the  same.  The  frequency  distribution  of  many 
such  locations  might  be  as  indicated  in  figure  2.2. 
Here,  the  horizontal  and  vertical  coordinates  are 
the  x  axis  and  frequency,  respectively.  The  un¬ 
known  but  true  location  of  the  1010-mb.  isobar 
is  at  x*.  Areas  under  the  curve  represent  prob¬ 
abilities;  the  total  area  is  unity.  Thus,  the  prob¬ 
ability  that  the  isobar  will  be  located  within  a 
specified  interval,  i  Ax,  of  its  true  location,  x*  ,  is 
represented  by  an  area  on  the  graph.  The  shad¬ 
ed  area  is  an  example. 


Figure  2.2.  Hypothetically  True  Frequency  Distribution  of 
Isobar  Location. 

Naturally,  in  practice  such  a  probability 
cannot  be  found  because  repeated  random  sam¬ 
plings  of  a  true  pattern  by  observation  networks 
are  not  available.  Yet,  an  upper  limit  to  this 
probability  can  be  estimated  from  available  in¬ 
formation.  This  is  accomplished  by  subjecting 
the  synoptic  analysis  itself,  rather  than  the  true 
pattern,  to  hypothetical  random  samplings  by 
networks  and  evaluating  the  results.  Details  of 
the  sampling  theory  will  not  be  given  here,  but 
the  method  for  computing  limiting  probabilities 
can  be  described. 


The  analysis  based  on  these  reports  places 


As  before,  assume  that  Ax.  represents  the 
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error  in  location  of  an  isobar  along  the  perpen¬ 
dicular  x  axis.  With  numerous  random  sam¬ 
plings  a  large  number  of  such  errors  is  obtained. 
Now,  assume  that  the  mean  of  the  squares  of 
these  errors  is  the  variance: 

*2(*M A?)2,  (1) 

where  the  bar,  here  and  hereafter,  represents 
the  arithmetic  mean  of  the  quantity  under  it. 
According  to  the  network  sampling  theory. 


How  are  these  results  of  use  in  synoptic  an¬ 
alysis?  A  numerical  example  will  illustrate. 
Suppose  that  a  sea-level  pressure  analysis  is 
made  over  an  area  of  6,250,000  square  miles, 
approximately  the  area  of  South  America.  There 
are  100  synoptic  pressure  reports  throughout 
the  region.  We  want  to  know  something  about 
the  accuracy  of  the  analyzed  pattern. 

With  R  -  6,250,000  square  miles  and  n  = 
100,  equation  (3)  yields 


o2(x)  =  0.056  l2, 


(2) 


l2  =  62,500  (miles)2, 


where  0.0  56  is  a  dimensionless  constant  and  i 
is  an  average  distance  between  synoptic  obser¬ 
vations.  The  latter  can  be  evaluated  from  the 
equation. 


l2  =  R /n,  (3) 

where  R  is  the  area  of  the  analysis  region,  and 
n  is  the  number  of  synoptic  reports  in  the  region. 
The  use  of  (3)  will  be  descussed  further  in  sec¬ 
tion  4.1. 


or 


l  =  250  miles. 


Then  from  (2), 


o2(x)  =  3,500  (miles)2, 


or 


o(x)  =  59.2  miles. 


Theoretically,  <72(x)  is  the  variance  of  a  nor¬ 
mal  or  Gaussian  frequency  distribution,  as  illus¬ 
trated  in  figure  2.3.  The  coordinates  are  the 
same  as  in  figure  2.2.  The  mean  of  the  distri¬ 
bution  is  x,  with  standard  deviation  a .  Again, 
areas  under  the  curve  represent  probabilities, 
and  the  total  area  is  unity.  The  probability  that 
deviations  from  x  will  not  be  greater  than  i  Ax 
is  represented  by  an  area  such  as  the  shaded  one 
in  this  figure. 


Let  Ax  =  100  miles.  From  the  ratio 


Ax  100 
^x)=593 


1.7. 


Refer  to  figure  2.4  which  is  a  graph  of  \x/o(x) 
versus  probability,  derived  from  tables  of  the 
normal  distribution  [1],  Enter  this  figure  with 
1 .7  and  find  that  the  ordinate  value  is  91  percent. 
This  is  the  theoretical  probability  that  at  any 
point  the  analyzed  position  of  an  isobar  lies  with¬ 
in  -  100  miles  of  its  hypothetically  true  position. 


Figure  2.3.  Model  Frequency  Distribution  of  Isobar  Loca¬ 
tion. 


Now,  compare  figures  2.2  and  2.3  which  are 
drawn  to  the  same  scale.  For  the  same  choice 
of  Ax  in  both  figures,  the  shaded  area  in  figure 
2.3  is  larger  than  the  one  in  figure  2.2,  meaning 
that  the  corresponding  probabilities  are  unequal 
in  the  same  sense.  When  Ax  =  100  miles,  the 
larger  probability  is  91  percent.  Theoretically, 
then,  the  probability  that  the  analyzed  position  of 
an  isobar  at  any  point  lies  within  100  miles  of 
its  actual  true  position,  x*  ,  is  less  than  91  per¬ 
cent.  We  do  not  know  where  x*  is,  but  we  do  have 
a  limiting-probability  statement  concerning  its 
proximity. 

In  this  example,  it  was  assumed  that  the 
probability  under  the  normal  curve  exceeded  that 
under  the  "true”  curve  (fig.  2.2),  for  the  same 
Ax.  This  is  a  general  prediction  of  the  network 
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_ OP 

o-(x)  o-(q) 


Figure  2  4.  Cumulative  Probability  Versus  Ratios  Ax/a(x)  and  Aq/o(q). 


sampling  theory;  the  limiting  probability  should 
be  at  least  as  large  as  the  true  one.  Whether 
this  happens  in  reality  can  only  be  determined 
by  tests.  Several  tests  have  been  run,  all  suc¬ 
cessful  so  far.  Examples  will  be  given  later. 


is  probably  somewhat  different. 

What  is  the  relationbetween pressure  error, 
Ap  ,  and  isobar-location  error.  Ax  ,  at  a  given 
point?  From  the  calculus. 


It  should  be  noted  that  the  choice  of  Ax  =  100 
miles  was  arbitrary.  Other  values  are  possible; 
for  example,  50  miles.  Then 


Ax  50 
a(x)  59.2 


0.845 , 


and  the  limiting  probability  is  60  percent,  as  can 
be  verified  from  figure  2.4. 


The  accurate  location  of  an  isobar  may  not 
be  of  great  interest,  but  equations  (1),  (2),  and 
(3),  and  figure  2.4  are  fundamental  to  develop¬ 
ments  of  more  interesting  and  advanced  topics 
discussed  below. 


2.2  Specification  of  Pressure 

Associated  with  errors  in  location  of  isobars 
are  errors  in  pressure.  For  example,  at  the 
arbitrary  point  A  in  figure  2.1  there  is  a  1010- 
mb.  analyzed  pressure,  whereas  the  true  value 


(4) 

where  dp/dx  is  simply  the  gradient  of  p  in  the  x 
direction.  Note  that  the  x  axis  was  chosen  along 
the  gradient. 

The  true  gradient  of  p  at  point  A  is  not  known, 
but  it  can  be  approximated  from  the  analysis. 
Our  procedure  now  is  similar  to  that  of  the  pre¬ 
ceding  section;  square  and  average  errors  to  ob¬ 
tain  variances.  Thus,  the  variance  of  p  becomes 

=  (5) 

In  this  operation,  the  value  of  dp/dx  is  to  be  re¬ 
garded  as  constant. 

As  before,  o2(x)  is  considered  to  be  the  var¬ 
iance  of  a  normal  distribution.  Then  because 
dp/dx  is  a  constant,  a2( p)  in  (5)  is  also  the  vari¬ 
ance  of  a  normal  distribution,  according  to  sta¬ 
tistical  theory.  Combination  of  (2)  and  (5)  yields 
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o2(p)  =  0.056  l2( .  <6> 

To  illustrate  the  practical  use  of  (6),  we 
employ  the  example  of  figure  2.1.  Assume  that 
the  pressure  gradient  measured  at  A  from  the 
analysis,  is  3  millibars  per  100  miles,  and  that 
the  average  distance  between  synoptic  reports 
again  is  250  miles.  Then,  from  (6), 

°2(p)  =  3.15, 
or 

o(p)  -  (3.15)*  -  1.77. 

Select  a  value  of  Ap  =  2  mb.  Form  the  ratio 
mLn- 

Enter  the  graph  in  figure  2.4  with  the  value  1.13 
and  find  the  corresponding  probability,  74  per¬ 
cent.  This  means:  the  probability ,  that  the  analyzed 
pressure  at  A  lies  within  ±2  mb.  of  the  true  pressure  at  that 
point,  is  not  greater  than  74  percent .  Under  the  cir¬ 
cumstance  that  synoptic  reports  are  widely  scat¬ 
tered,  a  stronger  statement  than  this  hardly 
seems  possible. 

2.3  Comporison  of  Theory  and  Observotion  in  Pressure 

Specification 

A  probability  statement  like  the  one  above 
cannot  be  tested  at  a  single  point  in  an  analysis. 
Observed  probabilities  are  needed  for  compar¬ 
ison  with  theory.  Because  actual  synoptic  sit¬ 
uations  do  not  repeat  themselves,  observedprob- 
abilities  at  a  fixed  point  are  not  available  and 
it  is  necessary  to  make  indirect  comparisons, 
one  of  which  will  be  described  here. 

A  synoptic  situation  was  selected  at  random: 
the  sea-level  pressure  for  the  entire  Northern 
Hemisphere,  1230  Greenwich  mean  time,  19 
February  1950.  There  were  718  available  sy¬ 
noptic  reports  of  pressure  and  many  of  these 
were  accompanied  by  observed  wind  values.  Of 
these  reports,  3  53  were  deleted  on  a  random 
basis  and  the  remaining  3  65  were  plotted  on  a 
base  map  and  analyzed  without  prior  knowledge 
of  the  isobaric  pattern.  Where  available,  report¬ 
ed  winds  aided  the  analysis.  Then  the  actual  er¬ 
ror  in  pressure  at  each  of  the  3  53  deleted  sta¬ 
tions  was  found  by  comparison  of  analyzed  and 
true  values.  The  theoretical  variance,  a2(p)  , 
was  computed  at  each  of  these  locations,  by  use 
of  (6).  (The  same  value  of  i,  634  miles,  appear¬ 
ed  in  all  computations.  This  value  was  found 
from  (3),  with  n  =  3  65  and  R  the  area  of  the 


Northern  Hemisphere.) 

Actual  errors  and  theoretical  variances 
were  combined  into  separate  frequency  distribu¬ 
tions,  but  the  details  of  this  step  are  omitted 
here.  Both  distributions  are  shown  in  figure  2.5. 
The  smooth  curve  represents  a  theoretical  nor¬ 
mal  distribution,  while  the  broken  curve  is  a 
histogram  of  observed  errors.  The  important 
feature  is  that  the  theoretical  curve  shows  less 
dispersion  from  the  mean  (set  at  zero)  than  does 
the  observed  curve.  In  other  words,  the  observ- 
edprobability  of  deviations  to  lie  within  any  spe¬ 
cified  interval  centered  about  the  mean  is  not 
greater  than  the  theoretical  probability  for  the 
same  interval.  This  conclusion  supports  the  net¬ 
work  sampling  theory,  and  is  typical  of  results 
obtained  from  other  similar  studies. 


Figure  2.5.  Observed  Frequency  Distribution  (histogram) 
and  Theoretical  Frequency  Distribution  (sym¬ 
metrical  curve)  of  Pressure- Analysis  Errors  for 
Northern  Hemisphere ,  1230  GMT,  19  February 
1950. 


2.4  Extropolotion 

In  previous  sections,  the  primary  concern 
is  the  probabilities  of  accurate  synoptic  analy¬ 
sis.  Henceforth,  attention  will  be  focussed  on 
probabilities  of  accurate  predictions  from  anal¬ 
yses.  We  begin  with  simple  extrapolation. 

Figure  2.6  shows  three  hypothetical  loca¬ 
tions  of  an  isobar  segment  along  a  perpendicu¬ 
lar  x  axis.  The  past,  present,  and  extrapolated 
positions  are  at  xQt  xlt  and  x2,  respectively;  with 
distances  between  positions  being  equal.  Thus, 


(7) 
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Because  of  analysis  errors,  the  displacement 
x  -  xr  is  somewhat  in  error,  and  thereby  causes 
an  error  in  the  future  displacement,  x2  -  Xj : 

A  (x 2  —  x  ^  ~  x0)  *  (8) 

This  can  be  rewritten  in  terms  of  the  loca¬ 
tion  errors  themselves: 

Ax2  -  Ax  j  —  Ax j  -  Ax0  ,t 
or 

Al2  =2  Al, -Al0.|.  (9) 

Here,  the  error  in  predicted  location,  Ax2, 
is  expressed  in  terms  of  analysis  errors.  Now 
both  sides  of  (9)  are  squared  and  averaged  to 
obtain  a2(x2 ),  the  variance  of  prediction  errors. 
Thus, 


a2( x2)  =  0.056  i2  (4  +  1)  =  0.280  l2  .  (11) 

Now,  suppose  t  =  2  50  miles  as  before.  From 

(ID, 

<t2(x2)  =  17,500  (miles)2 
or 

o(x2)  =  132.29  miles. 

Select  an  interval  Ax2  =  100  miles.  Form  the 
ratio 

Ax,  100 

a(x2)  =  132.29  =  0  76  ' 

With  figure  2.4,  this  yields  55  percent  prob¬ 
ability.  In  a  network  whose  average  distance  between 
stations  is  250  miles,  the  probability  that  an  isobar  extrap¬ 
olated  by  use  of  equation  (7)  will  lie  within  1 100  miles  of 
its  true  location  is  55  percent,  at  best. 

Two  features  of  this  example  should  be  no¬ 
ted.  It  has  been  implicitly  assumed  that  extrap¬ 
olation  is  a  perfect  prediction  method,  because 
no  method  errors  were  considered.  However, 
the  italicized  statement  above  still  holds,  even 
if  extrapolation  is  a  poor  procedure. 

Secondly,  the  results  are  not  restricted  to 
isobar  displacements.  Thus,  x2  could  represent 
the  extrapolated  position  of  a  front,  isotherm, 
or  other  synoptic  feature. 

2.5  Pressure  Forecast  Using  Pressure  Tendency 


<r2(x2)  *  (Alj)2  =  4  (Al,)J  +  (Al0)J  -  4  Ax,Ax0  .  (10) 

We  will  assume  that  present  and  past  synoptic 
analyses  were  made  independently  of  each  other, 
or  else  that  the  time  interval  between  successive 
maps  is  large  enough  so  that  errors  Ax0  and  Ax} 
are  essentially  independent  of  each  other.  This 
makes  it  possible  to  set  the  cross-product  term 
AXjAxq  =  0,  as  the  two  errors  are  uncorrelated. 
In  practice,  either  of  these  assumptions  is  not  a 
severe  restriction.  And  it  is  most  important  to 
eliminate  Ax,Ax0  from  (9),  otherwise  it  cannot  be 
evaluated.  See  section  4.2  for  further  discussion 
of  this  subject. 

Next,  (AXj)2!  and  (AxQ)2[  are  each  replaced 
with  0.0  56  £2’,  from  (1)  and  (2),  so  that  (10)  be¬ 
comes 


The  pressure  at  a  fixed  point  can  be  fore¬ 
cast  from  the  relation, 

Pt  =  P0  +  bot-  (12) 

where  tlis  time;  b0  the  initial  tendency,  dp/dt; 
and  p0  and  p{  are  initial  and  forecast  pressures, 
respectively.  In  figure  2.7,  let  A  be  such  a  point, 
and  let  x2  and  x^axes  be  drawn  perpendicular 
to  isobars  (solid  curves)  and  isallobars  (dashed 
curves)  through  the  point,  as  an  example.  The 
forecast  error,  Apt  ,  depends  on  analysis  errors, 
Ap0  and  Ab0 1.  Thus, 

Ap,  =  Ap0  +  tAb0  .  (13) 

From  (4), 


(14) 
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and,  by  similar  reasoning. 


dbn 


Abo=XAX2 


where  db0/dx2  is  the  isallobaric  gradient, 
stitution  of  (14)  and  (15)  into  (13)  yields 


(15) 

Sub- 


APt=^fAx, 


db„ 


(16) 


The  variance  of  forecast  errors,  &2(pt),  is 
found  by  squaring  and  averaging  both  sides  of 
(16): 


Here,  dp Q/dx x  and  dbjdx2  are  treated  as  con¬ 
stants  to  be  evaluated  from  the  synoptic  analy¬ 
sis,  and  the  last  term  is  set  equal  to  zero  be¬ 
cause  analysis  errors  of  isobars  and  isallobars 
are  mutually  independent  (Ax  t  Ax2  =  0). 

For  a  numerical  example  of  (17),  let  i  = 
250  miles,  dpjdxl  =  3  mb.  / 100  miles,  and  dbjdx 
=  1  mb./3  hr.  per  100  miles.  Then,  as  in  equa¬ 
tion  (11), 

aHft)  “  0-35  (9  +  0.11 12) , 


A 


db„ 


-2(p,Map,)2  =  (IT)  (A*,)*  ,+  **(■£*)  (&x2) 


dPo  dbo 


(17) 


+  2tdTIlf2AX  1AX2 


Figure  2.7.  Example  of  Isobar  (solid  lines)  and  Isallobar 
(dashed  lines)  Patterns. 


where  t  is  measured  in  hours. 

Let  Apt  =  2  mb.  By  use  of  (18)  and  figure 
2.4,  limiting  probabilities,  P  ,  for  pt  can  be 
found,  as  in  previous  sections.  Results  for  var¬ 
ious  values  of  t  are  given  in  table  2.1.  (The  last 
column  shows  values  of  pt  predicted  from  (12), 
with  assumed  initial  values,  P0  =  1012  mb.  and 
b0  =  +  2  mb./ 3  hr.,  at  point  A.  )  Limiting  prob¬ 
abilities  of  successful  prediction  of  pressure 
within  -  2  mb.,  decrease  with  increasing  fore¬ 
cast  periods,  which  is  consistent  with  experi¬ 
ence. 


Table  2.1.  Predicted  Pressures  (pt)  and  Limiting  Probabil¬ 
ities  (P)forPressures  to  Be  Within  2  mb.  of  True 
Values  at  Various  Times  (t). 


t 

(hour) 

2 

»(pt) 

P 

(percent) 

p« 

(mb.) 

0 

1.128 

74 

1012 

6 

0.936 

65 

1016 

12 

0.676 

50 

1020 

18 

0.504 

39 

1024 

24 

0.396 

31 

1028 
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3.  SUMMARY  OF  METHOD 


With  the  examples  of  the  previous  section 
as  an  introduction,  we  can  outline  the  method  for 
obtaining  limiting  probabilities  quite  generally 
and  formally.  Some  comments  about  the  exam¬ 
ple  in  section  2.5  are  appended  to  illustrate  each 
step: 

1.  Select  variable,  q,  to  be  studied. 


2.  According  to  (12),  pt  is  a  function  of  p0, 
b0,  and  t  : 

P,  =  P,(Po<  bo-  0  •  (19) 

3.  From  the  calculus  and  (19),  the  total  dif¬ 
ferential  is  written 


2.  Express  this  variable  as  a  function  of 
other  variables. 


dpt 

Apt=^Ap 


9P, 


A  b„ 


aP, 

dt 


At 


(20) 


Obtain  total  differential  of  the  variable 
in  terms  of  other  variables. 


The  partial  derivatives  in  (20)  can  be 
found  by  differentiation  of  (12).  Thus, 


4.  Simplify  total  differential  by  omitting 
from  it  any  difference  terms  not  due  to 
synoptic  analysis  errors. 

5.  Repeat  steps  2,  3,  and  4  for  each  of  the 
other  variables  in  the  total  differential. 
Continue  the  repetition  for  all  new  var¬ 
iables  introduced  during  the  process,  as 
far  as  possible. 

6.  Substitute  all  total  differentials  obtained 
in  step  5  into  first  total  differential.  This 
gives  the  basic  error  equation. 

7.  Square  and  average  both  sides  of  the  bas¬ 
ic  error  equation.  This  gives  the  variance 
equation. 

8.  Set  equal  to  zero  the  cross-product 
terms  that  are  averages  of  products  of 
uncorrelated  errors. 


9.  Replace  each  (Ax^2  with  0.056  t2  . 


*P0 


I 

"  Sb~  St 


(21) 


When  (21)  is  substituted  into  (20),  the 
latter  becomes 


Apt  =  Ap0  +tAb0tt  b0At  .  (22) 


4.  The  last  term  in  (22)  can  be  omitted  be¬ 
cause  there  is  no  error  in  t  introduced 
by  synoptic  analysis;  the  forecast  period 
can  be  specified  accurately,  indepen¬ 
dently  of  any  analysis.  Therefore.  At  = 
0  and  (22)  reduces  to  (13). 

5.  The  other  variables  in  (13)  are  p0  and 
bQ.  Their  total  differentials  are  given 
simply  by  (14)  and  (15).  The  new  vari¬ 
ables  introduced  by  (14)  and  (15)  are 
and  x2  .  But,  the  latter  are  independent 
variables,  so  the  repetition  of  steps  2,  3, 
and  4  terminates  here. 


10.  Evaluate  E2  from  R/n,  and  partial  deriv¬ 
atives  from  ratios  of  finite  differences, 
using  data  from  synoptic  analyses. 

11.  Use  computed  value  of  a(q)  and  specified 
value  of  A q  to  find  limiting  probability 
from  figure  2.4. 

Comments: 


6.  This  step  results  in  (16). 


7.  This  gives  (17). 


8. 


The  term, 
from  (17). 


dPodbo 
dxx  dx2 


is  dropped 


9  and  10.  These  steps  lead  to  (18). 


1.  In  the  example  of  section  2.5,  we  let  q  =  pt. 


11.  See  table  2.1. 
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4.  DISCUSSION 


4.1  Evoluotion  of  t 

According  to  (3),  the  square  of  the  average 
distance  between  synoptic  reports  is 

£2  =  R/n  , 

where  R  is  the  geographic  area  over  which  the 
analysis  is  made,  and  n  is  the  number  of  reports 
in  that  area.  In  all  studies  to  date,  it  has  been 
satisfactory  to  compute  £  in  this  way.  For  ex¬ 
ample,  the  successful  results  described  in  sec¬ 
tion  2,3  were  obtained  by  letting  R  be  the  total 
area  of  the  Northern  Hemisphere  and  n  the  num¬ 
ber  of  hemisphere  stations  used  in  analysis. 

A  special  problem  arises  when  analyses  are 
made  over  oceans  where  data  tend  to  be  sparse. 
If  data  are  too  sparse,  no  reliable  analyses  can 
be  made.  If  there  are  some  data  permitting  an¬ 
alyses  of  large-scale  features,  at  least,  then  R 
should  be  made  the  ocean  area  although  m  is  a 
small  number.  However,  if  the  weather  systems 
to  be  analyzed  have  emerged  from  a  nearby  land 
area  having  a  dense  observation  network,  so  re¬ 
cently  that  by  continuity  the  analyst  still  has  a 
good  concept  of  their  features,  then  it  is  suitable 
for  computation  purposes  to  use  a  value  of  t 
characteristic  of  the  dense  land  network.  It  is 
important  here  to  realize  that  the  chief  goal  is 
a  limiting  or  maximized  probability  value,  and 
that  this  is  not  unique  in  a  given  problem.  This 
is  because  all  probabilities,  from  unity  down  to 
computed  limiting  value,  are  themselves  limit¬ 
ing  probabilities.  Of  course,  it  is  desirable  to 
obtain  a  value  as  close  to  a  lower  “true”  prob¬ 
ability  as  possible.  But  at  least  it  can  be  ensured 
that  a  probability  will  be  a  limiting  one  if  £  is 
made  sufficiently  small.  Network  densities  on 
land  are  suitable  for  this  purpose. 


guarantee  independence  of  errors.  An  apparent¬ 
ly  safe  guide  in  practice  is  that  the  map  interval, 
t,  satisfies  the  inequality 

'*f/Cn,a,  •  (23) 


where  c  max  is  the  maximum  apparent  speed  that 
synoptic  features  (isobars,  fronts,  etc.)  are  ob¬ 
served  to  move.  Six  hourly  and  twelve  hourly 
maps  satisfy  this  criterion  in  dense  networks, 
particularly  on  land. 


The  problem  of  cross-product  terms  also 
is  present  when  evaluating  errors  of  gradients. 
For  example,  suppose  it  is  desired  to  study  er¬ 
rors  of  the  east-west  component  of  the  geo- 
strophic  wind,  u,  on  a  constant  pressure  surface. 
Now, 

“-"ft*  (24) 


where  g  is  the  acceleration  of  gravity,  f  the 
Coriolis  parameter,  and  dZ/dy  is  the  north -south 
height  gradient  of  the  constant  pressure  surface. 
In  practice,  the  gradient  would  be  computed  as 
a  ratio  of  finite  differences: 


dy  L 


(25) 


where  L  is  a  selected  distance  along  the  y  axis 
centered  at  the  point  where  u  is  to  be  found,  and 
Zj  and  Z2  are  heights  interpolated  from  the  con¬ 
tour  analysis  at  end  points  of  L.  Then  (24)  be¬ 
comes 


8  (Z,  -  Z,) 
f  L 


(26) 


Following  the  general  method  in  chapter  3,  we 
regard  u  as  a  function  of  g,  f,  Zlt  Z2,  and  L  and 
write  its  total  differential  as 


4.2  Treatment  of  Cross-Praduct  Terms 


This  subject  arose  earlier.  In  section  2.5, 
the  errors  of  isobaric  and  isallobaric  analyses 
were  regarded  as  independent  of  each  other, 
which  is  a  reasonable  assumption,  so  that  the 
term  containing  AxxAx2  in  (17)  could  be  set 
equal  to  zero. 

In  the  extrapolation  example  of  section  2.4, 
the  justification  for  dropping  AXjAx0  from  (10) 
is  that  the  analyses  at  different  map  times  are 
made  independently  of  each  other  or  that  the  in¬ 
terval  between  these  maps  is  large  enough  to 


a,,  dw  * _  du  Af  du  aT 

Au“-Sg  A8  +  5f  Af-+az,AZi 


dU  .  du  A  f 

hdZ  AZ2  +ST  AL 


(27) 


The  underlined  terms  become  zero  because  g 
is  a  known  constant,  L  is  predetermined  exactly, 
and  f  depends  on  latitude  and  is  assumed  con¬ 
stant.  By  partial  differentiation  of  (26),  (27)  re¬ 
duces  to 

gAZ,  g(AZ2-AZ,) 

L -  •  (29) 


g  &Z2 

AU--J  — 


f  L  m~f 


The  total  differentials  of  Zl  and  Z2  are 

az, 

AZ1  =1x7  AX1 


and 


(29) 
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when  h  is  large  enough  (as  in  this  case)  to  be  of 
practical  value  or  interest.  The  variance  equa¬ 
tion  becomes 


AZ 


£2 
2  *  dX 


2  > 


where  axes  xl  and  x2  are  defined  to  be  along  the 
height  gradients  at  the  points  where  and  Z%. 
are  interpolated  from  the  analysis.  By  combi¬ 
nation  of  (28)  and  (29),  we  have 


A  0  <*Z2  <*Z, 

=  AX,) 


fL  ^x2 


(30) 


The  variance  of  u  is  obtained  by  squaring  and 
averaging  (30): 

az, 


<72(u)=(Au)  =(- 


£Z, 

ax 


£)(ax2)  ♦  (- 


az,  az, - n 

2w,  ir2  Ai>Axd  • 


(31) 


The  last  term  is  a  cross  product  that  can  be  set 
equal  to  zero  if  Axl  and  Ax2  are  uncorrelated. 
The  latter  condition  will  be  met  if  the  distance 
chosen  for  L  is  sufficiently  large,  according  to 
the  inequality 


=  0.056  e2  [(|L)2  +  (|L.)2] 


(37) 


after  the  cross  product  term  is  dropped. 


4.3  Presentation  of  Results 

Each  of  the  previous  examples  is  concerned 
with  a  single  probability  value  at  an  initial  or 
forecast  time.  Depending  on  the  meteorological 
variable,  it  may  be  possible  to  obtain  a  field  of 
predicted  probabilities  associated  with  the  field 
of  the  predicted  variable.  Thus,  suppose  the 
height,  Z  ,  of  the  500 -mb.  surface  is  to  be  fore¬ 
cast  at  numerous  grid  points  by  a  standard  meth¬ 
od.  Then,  at  each  point  the  forecast  variance 
°2(Z)  can  be  obtained  and  used  to  find  a  limiting 
probability  there.  The  grid  field  then  can  be  an- 
alayzed  with  isolines  of  probability. 


L  ll.  (32) 

This  is  not  a  severe  restriction  in  practice.  Then 
with  the  aid  of  (1)  and  (2),  (31)  becomes 

a  2  f  <9Z ,  2  aZ„  2] 

o2(u)  =  0.056  e2  (j£)  [(^)  ♦(^)j  .  (33) 

which  can  be  evaluated  with  the  aid  of  (3). 

As  a  final  instance  wherein  cross  products 
appear  but  can  be  eliminated,  consider  the  prob¬ 
lem  of  errors  in  determination  of  thickness  be- 
tween  two  isobaric  surfaces.  Let  h  be  the  thick¬ 
ness  between  the  700-mb.  and  1000-mb.  surfaces 
at  a  given  geographic  location,  so  that 

h  =  Z  ~  z  ,  (34) 


Figures  4.1  to  4.3  provide  an  example.  The 
first  figure  shows  the  analyzed  height  field  at 
500  mb.,  0000  Greenwich  mean  time,  8  January 
1963.  By  means  of  a  barotropic  model  which 
conserves  absolute  vorticity,  a  numerical  24- 
hour  forecast  was  made  from  this  initial  state, 
and  is  shown  in  figure  4.2.  Then,  the  method  of 
chapter  3  was  applied  to  the  barotropic-model 
equation  and  the  synoptic  features  of  figure  4.1 
to  yield  numerical  predictions  of  o2(Z)  at  grid 
points  throughout  the  North  American  region. 
Corresponding,  limiting  probabilities  were  plot¬ 
ted  and  analyzed  in  figure  4. 3.  The  value  of  each 
isopleth  is  the  limiting  probability  for  predicted 
height  values  to  lie  within  -  100  feet  of  true 
values  at  0000  Greenwich  mean  time,  9  January 
1963.  The  map  shows  areas  where  forecasts 
have  good  and  poor  likelihoods  of  success. 


where  Z  and  z  are  the  heights  of  the  upper  and 
lower  surface,  respectively.  By  standard  pro¬ 
cedure, 

dZ  dz 

Ah  =  AZ  -  AZ  m—  AXj  ,  (35) 

where  the  Xj  and  x2  axes  lie  along  the  height 
gradients  at  700  and  1000  mb.  Then, 


At  this  point  the  reader  may  have  the  follow¬ 
ing  question  in  mind.  Suppose  two  analysts  in¬ 
dependently  drew  separate  contour  fields  for  the 
synoptic  situation  of  January  8,  Their  analyses 
would  be  slightly  different,  consequently  their 
analyzed  probability  fields  for  January  9  would 
differ  from  each  other  to  some  extent.  Would 
both  fields  still  give  valid  limiting  probabilities? 


a2  (h)  =  (Ah)2  =  (f-)2  (AX  ,)2  +  (pf  (Axf 

~2axlaT2  axiax2  ■ 


(36) 


Now,  if  the  contour  analysesJhave  been  made  in¬ 
dependently  we  can  set  AxtAx2  =  0,  particularly 


From  tests  so  far  and  from  deliberate  intent 
in  development  of  the  method,  the  answer  is  yes. 
Synoptic  analyses  are  gross  smoothings  of  re¬ 
ality  that  presumably  resemble  each  other  more 
closely  than  anyone  of  them  would  resemble  the 
true  situation  in  all  its  complexity. 
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Figure  4.1.  Initial  Height  Field  (100’s  of  ft.)  at  500  mb.,  0000  GMT,  8  January  1963. 


Figure  4.2.  Forecast  Height  Field (100#s  of  ft.)  at  500  mb.,  0000  GMT,  9  January  1963. 
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Figure  4.3.  Forecast  Probability  Field  (percent)  at  500  mb.,  0000  GMT,  9  January  1963. 


4.4  Verifications 

A  comparison  between  observed  and  theo¬ 
retical  probabilities  was  given  earlier  in  figure 
2,5.  Another  type  of  comparison  which  has  a 
simple  interpretation,  will  be  described  in  this 
section. 

Suppose  that  forecast  and  verification  height 
fields  at  the  500-mb.  level,  at  a  specified  time, 
are  to  be  compared  with  each  other.  At  each  of 
N  points,  the  error  in  the  forecast  height  can  be 
computed.  A  frequency  distribution  of  the  N 
errors  will  then  reveal  what  fraction  fell  with¬ 
in  *100  feet  of  verification  values,  for  example. 

Next,  the  forecast  variance,  o2(Z) ,  at  each 
of  the  N  points  can  be  used  to  find  the  limiting 
probability,  Pt  ,  for  the  predicted  height  to  fall 
within  *  100  feet  of  the  true  height.  Then  one 
can  determine  the  average  (expected)  probabil¬ 
ity,  P  for  this  result  at  the  N  points,  from 

P-1f(P„PJ,...,Pt....,PN)  .  08) 

If  N  is  large  (greater  than  100),  and  P  equals 
or  exceeds  the  observed  fraction  of  errors,  this 
may  be  accepted  as  support  for  the  theoretical 
probability. 


Figure  4.4  shows  an  actual  comparison  be¬ 
tween  theoretical  probabilities  (T)  and  observed 
relative  frequencies  (0)  for  various  error  in¬ 
tervals  (M)s  This  is  a  verification  graph  for 
10  synoptic  cases  in  which  Z  was  predicted  by 
the  numerical  barotropic  model  mentioned  in 
the  last  section.  For  each  value  of  AZ,  N  =1,645 


Figure  4.4.  Observed (0) and  Theoretical  (T)  Limiting  Prob¬ 
abilities  for  Forecast  Heights  to  Be  Within  In¬ 
terval  AZ  of  True  Heights  at  500  mb.,  Based  on 
1645  Numerical  Forecasts  in  10  Synoptic  Situa¬ 
tions. 
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pairs  of  values  were  used.  It  is  apparent  that 
the  theoretical  probabilities  are  truly  limiting 
values  because  they  bound  the  observed  values 
from  above. 

Parenthetically,  it  may  be  noted  that  the 
predicted  variances,  a2(Z),  can  also  be  used  to 
determine  limiting  probabilities  of  predicted 
ge o st rophic  wind  components.  Thus,  for  the  u 
component,  the  variance  a2(u)  results  from 


squaring  and  averaging  (28),  to  give 

o2(u)  =  (^7  =  (JL)2  [(AZ/  +  (AZ/ J 

=  (A)2[>(Z2)  +  a’(Zx)]  ,  (39) 

from  which  the  cross  product  term  has  been  suit  - 
ably  eliminated.  Graphs  similar  to  that  of  figure 
4.4  have  been  computed  for  u  and  v  components, 
and  with  similar  success. 
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5.  APPLICATION  TO  EXTRAPOLATION  FORMULAS 


Given  below  are  five  equations  for  extrapo¬ 
lating  synoptic  features,  and  corresponding  vari¬ 
ance  equations  for  predicted  displacements. 
Each  method  refers  to  figure  5.1  which  shows 
three  consecutive  past  positions  and  the  extrap¬ 
olated  position  of  an  isoline  (  A,  B,  C,  and  X, 
respectively  along  an  x  axis.  It  is  not  necessary 
that  this  axis  be  straight,  as  shown,  but  it  should 
be  perpendicular  to  the  isoline  at  all  four  inter¬ 
section  points. 


Figure  5.1.  Successive  Locations  of  Synoptic  Feature  Along 
x  Axis. 


5.1  Constant  Speed 

Xi-C  +  l/t  ,  (40) 

where 


is  the  measured  speed  from  B  to  C  in  map  in¬ 
terval  r  ,  and  t  is  time.  Then 

<r2(x,)  =  0.056  [l  +2i  +2(i)2]  (42) 

is  the  variance  equation  for  Xj  . 

5.2  Varioble  Speed 

X2  =  C  +  vt  +J  at2  ,  (43) 

where  r  is  given  by  (41),  and 

a  (44) 

is  the  measured  acceleration  involving  distances 


B-A  and  C-B  ,  and  there  is  the  same  time  inter¬ 
val,  t  ,  between  successive  maps.  The  variance 
equation  takes  the  form 

*2(x2)  =  0  056  l2  [l  +  2i  +  3  (~)2  +  3  (l)3 

+  1-S^)4]  •  (45) 

5.3  Percentuoi  Chonge  Method  far  t  =  r 

x3-  C  =  (C  -  B)2/(B  -  A)  .  (46) 

*2(x3)  -  0.056  e2  (1  +  4k  +  8k2  +  4k3  +  2k1)  ,  (47) 

where 

fe  =  (C  -  B)  /(B  -  A)  ,  (48) 

here  and  hereafter. 

5.4  Percentual  Change  Method  for  t  «  2r 

x,  ~x3  =(x3  -Cf/(C-  B)  ;  (49) 

o2(x4)  =  0.056  l2  (1  +  8k  +  32k2  +  56fe3 

+  5 2k<  +  24ks  +  8k6)  .  (50) 

5.5  Wasko's  Method 

*5*“  C  =  C  -  B  +  ^  _  g  [(C  -  B)  -  (B  -  A)J  ,  (51) 

under  the  condition  that  C  -  B  ^  B  -  A. 

- - p —  (2  +  4k  +  12k2  -  4k3 

+  5fe4)  .  (52) 

References  [5,  7,  and  9]  can  be  consulted 
for  discussions  of  equations  (40)  through  (52). 
The  reader  may  find  it  useful  practice  to  verify 
the  variance  equations  by  the  procedure  of  chap¬ 
ter  3. 

For  practical  purposes,  the  variance  equa¬ 
tions  have  been  graphed:  (42)  and  (45)  in  figure 
5.2;  (47),  (50),  and(52)  in  figure  5.3.  Each  figure 
is  entered  along  the  horizontal  axis  to  find  a  val¬ 
ue  of  the  ratio  a(x)/i  on  the  vertical  axis.  Then, 
with  a  value  of  i  from  (3),  o(x )  can  be  found. 

For  an  example,  consider  the  extrapolation 
case  in  section  2.4.  Actually,  this  is  a  special 
case  of  method  5.1  when  t  =  r  .  Therefore,  enter 
figure  5.2  at  the  point  (t/r)  =  1  and  find  from 
curve  xx  that  (a/i)  =0.53.  Let  i  =  250  miles. 


then  a  =  132.5  miles,  which  is  close  to  the  value 
132.29  miles,  found  in  section  2.4. 
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t  /T 

Figure  5.2.  Values  of  o/i  for  Predicted  Displacements  xl 
and  x2  as  Functions  of  t/r. 
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Figure  5.3.  Values  of  a/i  for  Predicted  Displacements  xv 
x  vandx^{dashedline)  as  Functions  of  (C-B)/ 
(B-A). 
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